Numerical implementation of isolated horizon boundary conditions 
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We study the numerical implementation of a set of boundary conditions derived from the isolated 
horizon formalism, and which characterize a black hole whose horizon is in quasi-equilibrium. More 
precisely, we enforce these geometrical prescriptions as inner boundary conditions on an excised 
sphere, in the numerical resolution of the Conformal Thin Sandwich equations. As main results, 
we firstly establish the consistency of including in the set of boundary conditions a constant surface 
gravity prescription, interpretable as a lapse boundary condition, and secondly we assess how the 
prescriptions presented recently by Dain et al. for guaranteeing the well-posedness of the Conformal 
Transverse Traceless equations with quasi-equilibrium horizon conditions extend to the Conformal 
Thin Sandwich elliptic system. As a consequence of the latter analysis, we discuss the freedom of 
prescribing the expansion associated with the ingoing null normal at the horizon. 

PACS numbers: 04.25.Dm, 04.20.Ex, 04.70.Bw 



I. INTRODUCTION 

The general study of spacetimes containing a black 
hole whose horizon is in quasi-equilibrium is of direct 
interest in astrophysics and numerical relativity. A par- 
ticularly important application is the determination of 
inner boundary conditions for the construction of astro- 
physically realistic excised initial data for binary black 
holes in quasi-circular orbits 

QEHSHEi Beyond the 
construction of initial data, the analysis of these bound- 
ary conditions also provides a helpful insight into the evo- 
lution problem. More precisely, a generalization of these 
conditions constitute an integrant part of the modeliza- 
tion of black hole horizons as world-tubes of marginally- 
trapped surfaces, according to the characterizations in 
the quasi-local horizon formalisms of trapped and dynam- 
ical horizons (see reviews In this sense, the quasi- 
equilibrium conditions here discussed offer a test-ground 
for the general dynamical case, in a better controlled sce- 
nario. Lessons acquired in the quasi-equilibrium case can 
then be applied to the evolution of an excised black hole 
in a constrained scheme like the one proposed in Ref . . 
In addition, quasi-equilibrium conditions themselves are 
of direct interest in slow evolution schemes, like the min- 
imal no-radiation approximation proposed in Ref. |lfj| . 

The isolated horizon formalism developed by Ashtekar 



et al. (see Ref. Q for a review) provides a particu- 
larly well-suited framework for studying a black hole in 
equilibrium inside a generically dynamical spacctimc 1 . 
It is an example of trapped horizon in which the horizon 
world-tube is a null hypersurface. This null-like character 
encodes the quasi-equilibrium characterization. 

First derivations of quasi-equilibrium horizon bound- 
ary conditions were presented in Refs. 0,0 (see also Ref. 
Q). The detailed analysis of isolated horizons in a 3+1 
description of spacetime permits the systematization and 
extension of these results. Following this line of research, 
an ensemble of boundary conditions for the 3+1 fields 
has been proposed in Refs. 0, 0, In this paper 
we address the problem of testing numerically these sets 
of boundary conditions. Translating Einstein equations 
into a specific system of partial differential equations, 
for which the isolated horizon prescriptions become ac- 
tual analytical boundary conditions, requires the choice 
of a particular resolution scheme for the geometrical field 
equations. For concreteness (and also motivated by the 
fully-constrained scheme in Ref. [9| , in which constraints 
are solved at each time step), we focus here on the con- 
struction of initial data in a Conformal Thin Sandwich 
(CTS) approach E^. Making use of an excision tech- 
nique, i.e. removing a sphere S of coordinate radius rn 
from the initial spatial slice S, and imposing S to stand 
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1 Throughout the paper we abuse the language and use the ex- 
pression black hole spacetime to refer to a spacetime containing 



an isolated horizon inside, 
event horizon. 



without any mention to the notion of 
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as a space-like slice of an isolated horizon, an ensemble 
of inner boundary conditions for the CTS elliptic system 
is determined. Note that throughout the paper we use 
dimensionless physical and geometrical quantities which 
result from the corresponding dimensional quantities by 
rescaling with the appropriate power of r# . 

Regarding index notation, greek letters denote 
lorentzian indices and are mainly used in section^ where 
null and space-like geometries appear in the same con- 
text. Latin indices refer specifically to objects living on 
space-like slices and are used from section ITTT1 on. 



II. BOUNDARY CONDITIONS: 
GEOMETRICAL FORM 

A. Geometrical boundary conditions 

In order to introduce a characterization of quasi- 
equilibrium, we need a notion of time evolution. A nat- 
ural evolution vector on a null world-tube H, sliced by 
a given family of trapped surfaces {St}, is provided by 
the null- vector £ that Lie-draggs S t onto S t +st- We con- 
sider the horizon slicing {S t } as induced by a 3+1 space- 
time foliation of space- like surfaces Et. Firstly, we fix 
the notation. We denote by n a the time-like unit nor- 
mal to St, by N the associated lapse, and by s a the unit 
space-like normal to St lying on E t . The horizon evo- 
lution vector is then expressed as £ a = N ■ (n a + s a ). 
The ingoing null vector k a (in the plane defined by n a 
and s a and normalized as k^£^ = —1) is written as 
k a = ^(n a - s a ). We denote by (j a0 ,K al3 ) the 3- 
metric on S t and the extrinsic curvature (with sign con- 
vention K a p = -\Cnlu$ = ~Y a V >/?). The induced 
metric on the marginally trapped surface S is then given 
by q a = la/3 - s a sp. 

Quasi-equilibrium boundary conditions follow from pre- 
scribing certain 3+1 fields to be time independent on the 
horizon. In addition, other relevant boundary conditions 
(not necessarily related to quasi-equilibrium) follow from: 
i) analytical requirements on the well-posedness of the el- 
liptic system, ii) numerical control of the horizon slicing 
taking into account the geometry of the horizon and, iii) 
choice of coordinate system adapted to the horizon (ul- 
timately motivated by numerical reasons). We briefly 
review this ensemble of conditions (for their systematic 
derivations and justifications, see Refs. [H|ll (13). 



B. Quasi-equilibrium conditions 

Prescribing the time-independence of a particular com- 
bination of 3+1 fields can represent either an actual re- 
striction on the geometry of Ti. as a spacetime hypersur- 
face, or rather it can refer to the manner H, is described 
in the 3+1 slicing. Both cases are relevant in a numerical 
relativity context. 



1. The minimal notion of quasi-equilibrium is pro- 
vided by the non-expanding horizon condition, namely 
the time independence of the induced metric on S: 

1 ojt ' B^i^iMv = 0. Expressed in terms of the expansion 
9(1) and shear O(t) associated with the outgoing null nor- 
mal £ a , i.e. 

0( t) = eTV^ , (i) 

(o-(i))a/3 = q^ a q" p^^u - )^{i)q a p , (2) 

this amounts to fP. [TlLIT^ : see also Ref. for a heuris- 
tic discussion on the vanishing of the shear) 

9 w = = (a {i) ) afj . (3) 

These three conditions [a(£) is a symmetric trace-less ten- 
sor on S 2 } constitute an actual restriction on the geome- 
try of H, essentially related to its null-character via the 
Raychaudhuri equation. Physically, they mean that the 
area of the horizon remains constant in time. 

2. Another proposed quasi-equilibrium condition pl| 
consists in prescribing the time-independence of the ver- 
tical (i.e. in the £ a direction) component of the angular 
variation of the null normal £ a . Explicitly this translates 
into C(fl a = 0, where fi Q := — q fl a (fc„V ^l"). Recasting 
VL a in terms of 3+1 fields, this condition means that the 
combination fl a = 2 D a h\N — K^s^q" a remains constant 
in time (where 2 D is the connection associated with q a p). 
Its underlying justification uses the weakly isolated hori- 
zon notion (see Ref. |ll|). However, we can heuristically 
motivate it in two manners: 

i) Given St with an axial symmetry generated by </>^ and 
with volume element 2 e = ^Jqdx a f\dx b , an angular mo- 
mentum Ju = — g^jj J s 2 e can be associated with 
the horizon |17| . The surface density of this momentum, 
Q a , satisfies a Navier-Stokes-like evolution equation for a 
viscous fluid (see [l3l[T8l | for a discussion of each term): 

c e n a + en a = %^ a T^t + 2 D a n 

- 2 D^ a + l - 2 D a 9 , (4) 

where K is the non-affinity coefficient of the null geodesic 
generated by t. V ' tl a — k £ a . Non-expanding con- 
ditions in Eq. (0 imply the vanishing of the viscous 
terms in Eq. |0J (as well as the force surface density 
term; see [l3|)) which reduces then to the Euler equation, 
C(P. a = 2 D a K. A natural quasi-equilibrium condition in 
this fluid analogy is given by CiQ a — 0, which translates 
into: 2 D a n = 0. This expresses the constancy of k on 
S. Writing then k — k — const, its 3+1 decomposition 
provides an evolution equation for the lapse: 

dlnN = k - s^D^N + NK fiV s fi s" . (5) 

Under the gauge choice CiN = 0, this equation becomes 
the boundary condition proposed in Ref. |ll| : 

K = s' i D l ,N-NK^s v . (6) 
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This condition does not restrict the geometry (and there- 
fore the physical features) of the horizon, but rather 
chooses a convenient 3+1 description. 

ii) From the relation above CeQ a = 0, one can motivate 
the time-independence of £l a from the Hamiltonian anal- 
ysis of isolated horizons (see 0, El)- I n this context, 
the non-affinity coefficient k is interpreted as the surface 
gravity of the horizon. The constancy of k turns out to 
be a very natural equilibrium condition as the quasi-local 
zcroth law of black hole mechanics. Such a Hamiltonian 
analysis provides a canonical constant value for the sur- 
face gravity, namely k q — K Korr (a, J), where a and J are 
the area and angular momentum of St, and K Kcrr (a, J) is 
the corresponding surface gravity of a Kerr black hole. 



C. Other geometrical boundary conditions 

Together with quasi-equilibrium motivations, some ge- 
ometrical inner boundary conditions follow from genuine 
numerical motivations. As a first instance of such bound- 
ary conditions, in some numerical schemes it is important 
to keep the location of the horizon fixed at a given po- 
sition. Geometrically this means that the 3+1 evolution 
vector t a = Nn a + f3 a (where (3 a is the shift vector) 
must be tangent to the horizon hypersurface TL. De- 
composing the shift in its normal and tangential parts 
to S, i.e. f3 a = bs a - V a (with F% = 0), it follows 
t a = i a - V a + (6 - N)s a . Therefore t a is tangent to TL 
iff 



b- N = 



(7) 



A second example of numerically motivated boundary 
condition follows from the need of a prescription for the 
lapse that incorporates information on the geometrical 
content of the horizon, but still leaves some rescaling free- 
dom to control the magnitude of N. In other words, a 
geometrical prescription for the slicing {St} of TL, rather 
than for N itself. In Ref. |2fj it is shown that specifying 
2 DT2 A ( ) i.e. the divergence of the 1-form f2 M introduced 
above, fixes the slicing of TL in an intrinsic manner (see 
also Ref. 01 f° r a 3+1 discussion). From the 3+1 ex- 
pression of fl a it follows 



III. CONFORMAL THIN SANDWICH 
DECOMPOSITION 

A. CTS Equations 

Fixing a representative 2 7^ in the conformal class of 
jij, we perform the following decomposition of the 3- 
metric and the extrinsic curvature 



Hj 
where 

with -y" 1 := £tj v , and 



1 



K 13 = A 13 + -K 1 



1 

2N 



(LP) 13 + 7 



~>:i 



K = - ■■' K, 



(9) 



(10) 



(11) 



and Di the Levi-Civita connection associated with jij. 
Inserting this decomposition in the Hamiltonian and 
momentum constraints, and prescribing on E the time 
derivative of the trace of the extrinsic curvature (K := 
CtK), results in the CTS equations. This is an el- 
liptic system for ("J, (3 l , N), once the initial free data 

{^ij , K, 7 J , K) are prescribed on the initial S. If a sphere 
iS is excised in S, as it is the case here, the boundary con- 
ditions for P l , N) on S are also a part of the free data. 
For concreteness, in this work we focus on free data 3 

with 7 J = = K, and consider different possibilities 
for (^/ij,K). Vacuum CTS equations, after a conformal 
rescaling of the lapse TV = ^ a N, have the form 



Kj3 l + -D i D k /3 k + Rlf3 k - N- 1 (L/3) ik D k N 
3 

-(a - 6)%- 1 (Ll3) ik D k i$ = -^ a Nb z K 



0, 



(12) 



AN 
+N 



•2(o- 



a ~ 



- l)D k ln^D k \nN 

a- A 



32 



12 

4-2a aV-1 



1 + a(a + l)Z? fe ln1'L> fe ln* 



N- l (Lf3) tJ (L(3)^ = ^- a (3 k D k K 



Different choices of the rescaling exponent a have been 
considered in the literature: a = 6 for defining the con- 
formal lapse in Refs. [TiL ITr| . a = — 1 is used in the 



2 AlnN = 2 Df{q^ p K flv s u )+ 2 D^ tl , (8) 

which fixes IniV up to a constant. In consequence, the 
lapse is specified modulo a multiplicative constant that 
can be chosen to keep the slicing under numerical control. 



2 We use latin indices from now on to emphasize that equations 
are defined on a space-like slice S. 

3 This choice is used in the literature QJ 0. as a quasi- 
cquilibrium prescription in the bulk. Gauge horizon condition Q 
gains then geometrical meaning, since it links bulk and horizon 
quasi-equilibrium notions by making the evolution vector t a to 
be tangent to the horizon. 
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numerical implementations of Refs. U,|5j,|6j, and a = — 2 
in Ref. The problem we address in this work is the 
numerical study of the elliptic system (|12(l . when com- 
pleted with the inner boundary conditions in Sec. [HJ 



B. Inner boundary conditions 

Inserting the conformal decompositions 4 10 in the ge- 
ometrical boundary conditions @, one finds 

45* A* + A* 4 * = - ^J-(Lf3) il s i s j (13) 
2N 

+ (l - pij&P 

from the vanishing expansion 0, 0> HH EH 113 , and 



(14) 



= , 



from the vanishing of the shear 

dEl, where P* = - 
V 1 , Hij = q , A5j is the (conformal) extrinsic curvature 

of S as a hypersurface in £ (with trace i? = q kl Hki), 

and we have used 7 J =0. In case of enforcing boundary 
condition J7J), the second bracket cancels, and the van- 
ishing of the first line of Eq. (|14|) characterizes V % as a 
conformal Killing vector of (St, qij) pLlllj. This provides 
a Dirichlet condition for V % once the conformal isometry 
has been chosen. 

In addition to the geometrical or numerical motiva- 
tions for the boundary conditions, we must also consider 
at this point the analytical well-posedness of the elliptic 
problem. A procedure to establish the uniqueness of the 
solution of an elliptic equation (even though not straight- 
forwardly generalizable to an elliptic system) consists in 
making use of a maximum principle 22, 23] . This in- 
volves the control of the convexity of the functions we are 
solving for (in our case depending on the choice of the 
exponent a), and in particular on the signs of their ra- 
dial derivatives at the boundaries. From condition (|13ll . 
control of PD^ on the horizon demands the control on 
the sign and size of (L0)ijS t s :l , something problematic 
if only the Dirichlet condition is imposed. An al- 
ternative condition proposed in 12], in the context of 
a Conformal Transverse Traceless (CTT) decomposition, 
consists in prescribing 5 * 6 • K; L jS % s^ to satisfy 



H ■ K ln s l s 3 < 



(15) 



We also introduce the rescaled induced metric fjij on St'- q~ij 



iSi, with s l = * 2 s\ 



Notation here differs from that in Ref. IT3 . Tilded objects 
there represents physical quantities, whereas here they refer to 
conformal counterparts according to 151 . 



Denoting f± = \P 6 • KijS l s : > this is enforced as 
2s k D k b-bH = 3N^^ 6 fi- 2 D k V k 



2V le D s s 



NK 



(16) 



If condition (JJJ is not imposed, the vanishing of the shear 
must be fulfilled by choosing, in addition to V 1 as a con- 
formal Killing symmetry, free data such that the traceless 
part of Hij vanishes (umbilical condition). More gener- 
ally one could solve condition (|14f> as an equation for 
V 1 [l3|. Extending the well-posedness analysis from the 
CTT to the CTS case is an important issue. Difficulties 
are two- fold. On the one hand, as pointed out in [24| . 
signs in Eqs. (JT^J for the (Lf3)ij (L(3) lj terms in the * 
and N are problematic for applying a maximum principle 
argument. No obvious choice of a cures the problem. On 
the other hand, the inclusion of new non-linear coupled 
boundary conditions for N makes the analytical problem 
even harder. However, a strong motivation for bound- 
ary condition (|15fl follows from its close relation with the 
characterization of St as a future trapped surface, i.e. 
with 6fu\ < 0, via the identity 



Kij S l s> - K 



2N 



N6 



(k) 



NO 



(fc) 



(17) 



where k a = ^(n a — s a ). More importantly, condition 
(|15fl is not specifically tied to quasi-equilibrium; in fact, 
0(k) < is part in the very definition of dynamical hori- 
zons, in order to garantee the horizon area increase law 
0- Moreover, such a condition on the sign of 6qa per- 
mits to exclude certain pathologies in the evolution of the 
horizon, e.g. the appearance of self-intersections [25J of 
the surface St, something to be avoided during the non- 
merger phase of the black hole evolution. The present 
quasi-equilibrium context offers a controlled test-bed for 
studying this condition. Numerical techniques seem to be 
an appropriate tool for a first analysis of this problem. 



IV. NUMERICAL RESULTS 

We make use of two independent codes for solving Eqs. 
ifT2)l. both using pseudo-spectral methods: the first one 
employs the elliptic solvers described in Refs. I26II27I l2Sj 
and implemented in the C++ library LORENE|29|. The 
second code has been specifically designed for the pur- 
pose of this paper. It uses a single domain technique for 
solving elliptic boundary value problems in the exterior 
of an excised spherical shell. 

In order to determine the elliptic system, we complete 
Eqs. (|12|l with a specific combination of five of the bound- 
ary conditions reviewed in the previous section. Differ- 
ent possibilities arise, but all of them must incorporate 
Eq. I|13|) and Eq. i|14|) . For concreteness, in this paper 
we will restrict ourselves to axially symmetric excised 
spheres (with azimuthal symmetry <jf), and impose as a 
Dirichlet boundary condition for the tangential part of 
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V 1 


b = P k ~su 


N 


0(e) = 


v = n„ ■ (f>* 


b = N 
* 6 ■ Kij S l s j = fi 


k = const 
2 D k Q k = f a 



TABLE I: Methodological assignment of boundary conditions 
to constrained fields. We keep fixed conditions for *]/ and V 1 
and study different combinations of the conditions for b and 
N. 



the shift: V 1 — il ■ 4> l , with Vl a constant [vanishing of 
the shear requires then appropriate additional boundary 
conditions or suitable free initial data in order to can- 
cel the second line in Eq. {HJ]. Eqs. (TSJ and (EJ 
fix three of the five boundary conditions. The remaining 
two will be chosen among Eqs. ©, ©, © an d ©3>- 
At this stage of the analysis it is methodologically useful 
to interpret each condition as associated with a specific 
equation of the system, even if this makes no strict sense 
due to the (non-linear) coupled character of the boundary 
conditions. Tabled summarizes this strategy, followed in 
Refs. [MS El ^2- The numerical implementations show 
that, keeping fixed the conditions for \I/ and V 1 , the dif- 
ferent combinations of conditions for b and N actually 
lead to the existence of solutions, at least for appropriate 
ranges of the free parameters (and independently of 7^ 
and K). The solutions are generically unique or finite in 
number, and only very particular (non-generic) choices 
of the functions f\ and fi lead to the appearance of an 
infinite number of solutions. 



A. Combination (6 = N, k = const) 

This set of boundary conditions implements the van- 
ishing of the shear for any choice of 7^ , as well as the 
constancy of the surface gravity. In Ref. 11] it was pro- 
posed the Hamiltonian canonical choice k = « Ke „ (a, J) , 
i.e. 



s l D t N - NK i:j s l S' 



(18) 



where K Korr (a, J) is a functional of , j3 l , N) . How- 
ever, the resulting elliptic system presents some prob- 
lematic features. We consider firstly the spherically 
symmetric case, fi = 0. The elliptic system admits 
then an infinite number of solutions, since maximal slic- 
ings of Schwarzschild provide the 2-parameter Estabrook- 
Wahlquist family of solutions j^D] , as pointed out in Ref. 
Q. Fixing the coordinate radius of the excised sphere 
only fixes one of these parameters, leading to a degen- 
erated problem. The very nature of this degeneracy 
suggests the way out, since setting the surface gravity 
to a given constant n fixes the representative of the 
Estabrook-Wahlquist family. As a code test, we have 
confirmed numerically that the system (b = N, k = K ) 
determines a unique solution, leading to a well-posed 
problem. In addition, the a posteriori evaluation of the 



quantity n Kcrr (a, J) on the constructed solution, results in 
K Korr (a, J) = k . This is in agreement with the degener- 
acy of the system [b = N,k = K Korr (a, J)], meaning that 
the operators on the left- and right-hand sides of Eq. i|18|) 
become identical on the space of solutions of this elliptic 
system. We conclude that, in the spherically symmetric 
case, a well-posed problem is defined by imposing the co- 
ordinate system to be adapted to the horizon (b = N) 
together with prescribing a given constant value for the 
surface gravity. In this system, the Hamiltonian canon- 
ical value for the surface gravity k = K Korr (a, J) cannot 
be imposed as a boundary condition, but it is actually 
recovered in the solution. Regarding the range of pos- 
sible values of k , this parameter is bounded by below, 
k > 1/8. In particular, this lower bound is associated 
with the vanishing of the lapse function on the horizon. 
No maximum value exists for n . However, the quadratic 
growth of the lapse when k increases, makes N to reach 
very large values very rapidly. For this reason the con- 
vergence of our codes is limited to a maximum value of 

K . 

More generally, we have also found that the system 
(b = N, K = K ) is well-posed in the rotating case, f2 7^ 0. 
On the other hand, the situation regarding the canoni- 
cal value of the surface gravity changes with respect to 
the spherically symmetric case: the a posteriori evalua- 
tion of K Kt}ll (a,J) does not provide identically k . This 
means that, in case of existing a solution to the system 
[b = N, k = K Korr (a, J)], the problem is not infinitely de- 
generated. In particular such a solution only exists if, 
when screening the solutions obtained by prescribing k, 
there is a value k for which K Kcrr (a, J) = k holds. Our 
numerical implementations seem to rule out this possi- 
bility, as illustrated by Fig. ^ This figure (implemented 
in conformal flatness, for concreteness) shows the rela- 
tive difference between the prescribed value of k , and 
the evaluation of K Korr (a, J) in the constructed solution. 
As in the spherically symmetric case, a minimum value 
for k is found, whereas the growth of the lapse limits 
the upper values we can numerically implement. In this 
range, and for f2 Q = 0.06, a non-vanishing minimum dif- 
ference between n Q and K Korr (a, J) is actually found. A 
similar behaviour is found for bigger values of f2 , even 
though the numerical limitations prevent us from deter- 
mining the minimum. 

Therefore in the rotating case fl 7^ the system 
[b = N, k = « Korr (a, J)] seem also to be ill-posed, but 
for the non-existence of solutions rather than because of 
the presence of infinitely many. 



B. Extending the well-posedness analysis from 
CTT to CTS: * 6 ■ K ij s i s j versus (£) prescription 

One of the main goals of this work is the study of con- 
dition ©>J), (O in the CTS system. In Ref. 12] it was 
shown that this condition, together with (|13|l and l|14|l . 
defined a well-posed problem in the CTT construction of 
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FIG. 1: Relative difference between « Kerr (o, J) and 

7 as a 

function of K a , for solutions to (b = N, k = k ) for different 
values of Q in the rotating case. 



initial data. A natural question is to study how this re- 
sult extends when the elliptic system is enlarged with an 
additional equation for the lapse. In a strict sense, this 
question is not properly formulated, since its answer will 
depend on the fifth boundary condition "for the lapse" . 
The aim here is rather to assess if a qualitative conclu- 
sion (independent of the details of the fifth boundary 
condition), can be formulated about the possible range 
of values of \& 6 • KijS l s^ with a focus on the negative 
ones. Our intention is not to perform an academical (nu- 
merical) extension of the CTT analytical result: we are 
ultimately motivated by probing some technical issues 
that will arise in the dynamical regime of the horizon. 
Indeed, given the relation with 9,^ via Eq. (|17|l . and 
the need to control the sign of 0,^ in a dynamical hori- 
zon, it is fundamental to know the range of values we can 
actually prescribe on St- 

a) Non-rotating case: Cl = 0. Let us firstly consider the 
spherically symmetric case. In addition to the boundary 
conditions in the first two columns of Table [0 and in or- 
der to probe the possible values of ^ 6 • K i jS l s : > , let us 
prescribe -K^s-'V = A < 0, with A a negative constant 
on S. We then complete the elliptic system with differ- 
ent boundary conditions "for the lapse" , and for each of 
them we screen the (negative) values of A for which we 
can construct a solution. Finally for each found solution, 
i.e. for each value of A, we plot the dimensionless quan- 
tity (* 6 • KijS^^/iH) against (9 {fc) )/(H). In Fig. I^we 
present the resulting curve where, for completitude, an 
extension to positive values of A has been included. The 
curve proves to be independent of the "lapse" bound- 
ary condition and, even though numerically we can only 
reach finite values of A, in this spherically symmetric case 
an analytical expression can in fact be obtained (see Ap- 
pendix^). For this reason, a compactification with the 
function arctan has been implemented for plotting Fig. 



121 

The existence in Fig. |2j of a minimum S m in for 
\l/ 6 • KijS L s^ and of an asymptotic negative non- vanishing 
value Sasym < for fy 6 ■ KijS l s : > (when KijS % s^ — > — oo), 
show that: i) the possible negative values of 'J 6 • KijS l s^ 
are bounded by below: S m i n < ('J 6 • K^s 1 s J ) / '(H); and ii) 
there is a range 5 asym < ■ K^s^s^)/ (H) < for which 
the prescription of the value of ■ KijS l s^ determines a 
unique solution for the elliptic system, exactly as it was 
concluded analytically in the CTT case. Let us also note 
that uniqueness is lost when the value \I/ 6 • i^s's 3 is pre- 
scribed between the values 5 m in < ^ 6 ■ KijS l s^ < S aS ym, 
for which two solutions exist. This is just another exam- 
ple of the non-uniqueness issue associated with solutions 
of the CTS equations, already pointed out in Ref. |24j . 
This is in contrast with the prescription of a negative 
non-conformally rescaled KijS l s J = 9^y. in the spheri- 
cally symmetric case, for each negative value of 9, fa (not 
bounded by below) there exists a unique solution. This 
suggests 9,^ as the good function to be prescribed on 
S, since no knowledge additional to the sign of 9^ is 
needed for consistency (see below). 

The situation is reversed for positive values of A, where 
only small values for 9,^ can be prescribed, leading al- 
ways to a degenerate solution. The good parameter is 
then *5> 6 ■ KijS l s 3 . However, these solutions do not de- 
scribe a future marginally trapped surface and their in- 
terest in the present context is only formal. 

b) Rotating case: fl 0. The independence of the curve 
"f 6 ■ KijS % si vs. 9,u with respect to the "lapse bound- 
ary condition" disappears in the rotating case. This is 
illustrated in Figs. 03 El and , where different bound- 




-1.4 -1.2 -1. -0.8 -0.6 -0.4 -0.2 0.0 



arctan [^K^s's^ / Hj 

FIG. 2: Numerical analysis of condition 1)15^ . Values of 
= JCysV /H and * 6 ■ /\.. r - .- // in the spherically 
symmetric case, after compactification by means of the arctan 
function. 
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ary conditions for N have been implemented 6 (due to 
the angular dependence, we plot now the minimum of 
• /\,,.s'.s-';. 

The most dramatic qualitative change in the rotating 
case is the existence of a certain value of il such that, for 
bigger values of this rotation parameter, the curve does 
not pass through the origin (in the case of Figs. El and 
this is true for any Q Q ^ 0). Arbitrarily negative values of 
"f 6 • KijS l s : > are found for small values of 9,j,y Therefore 
in this range of and in contrast with the spherically 
symmetric case, the negative values of \1/ 6 • KijS l s^ are 
only bounded from above. This shows that conclusions 
in Ref. [H for the CTT case do not extend to the CTS 
one for arbitrary values of V 1 = fin • • Given the variety 
of behaviours in Figs. OHElandEl it is difficult to extract 
generic conclusions, i.e. independent of the boundary 
conditions on N, about the possible values of 9,^-. and 

'J 6 • KijS l s ; ! . Our numerical simulations do not provide 
a complete understanding, but rather some restricted in- 
sight, of the boundary conditions properties in the case 
of a rotating black hole. 

Having stated this clearly, we highlight Fig. [3] corre- 
sponding to the interpretation of boundary condition (J7|) 
as the "lapse" boundary condition, N = b. In this case, 
we can actually formulate a concrete statement about the 
good parameter to be prescribed: for all values of fl 0l 
there exists a (small) negative value KijS l s' J — Aji < 
such that the prescription of 9,^ to A < A fio deter- 
mines a unique solution. It is very difficult to determine 
numerically if Aq o /H is actually zero or a very small 
value depending on fi . At this value KijS % s^ — Aji the 
conformal factor seems to have a pole and the quantity 
'J 6 • KijS l si diverges to negative values. In consequence 
curves in Fig. [3] for f2 ^ do not reach the origin, 
as pointed out above, changing the qualitative behaviour 
discussed in the spherically symmetric case. Still there 
exists a critical Sl C rit{~ 0.15) such that, in the range 
Q a < £l C rit, all curves present a local minimum S m i n and 
an asymptotic value S asy m at KijS l s^ — > -co. However, 
uniqueness (if existence at all) is lost when prescribing 
'J 6 • K i jS l s : > > Sasym- because of the negative divergence 
of \E' 6 - K ijS z s : > , either there are two or none solutions. The 
ultimate reason for focusing on this "fifth" boundary con- 
dition, is that it presents some geometrical/physical ad- 
vantage with respect to the other ones: i) it implements 
the NEH condition for any choice of -yy without the need 
of relying on the umbilical condition (see foot note El and 
ii) by enforcing the evolution vector t a to be tangent to 
H, which is in quasi-equilibrium, it chooses a coordinate 
system on the horizon in which no time-dependence is 
artificially introduced as a gauge effect (see also footnote 
13 for actual physical consequences in the binary case). 



N = b boundary condition 



o - 



'Hi 



-3.5 




-3 -2.5 -2 -1.5 -1 -0.5 

FIG. 3: Values of * 6 ■ K i jS i s i /H and 8 {k) /H, for 
boundary conditions KijS 1 s 3 — A, b — N and Q = 
0, 0.10, 0.11, 0.12, 0.20. Curves depart from the spherically 
symmetric case (dotted curve corresponds to Q = 0) as the 
rotating parameter increases. Solutions exist for every pre- 
cribed 9,^ and are unique. 



ify-sV = A, N = 0.2 




-3 -2.5 -2 -1.5 -1 -0.5 



FIG. 4: Values of * 6 • K ij s i s j /H and 6 (k) /H, for 
boundary conditions KijS 1 s 3 = A, N = 0.2 and Q. = 
0.00, 0.01, 0.02, 0.15. If the rotating parameter is smaller 
than Sl w 0.12620, the curves reach the coordinate origin, 
otherwise they diverge to — oo as 9,f,-./H — » 0. 



In sum, we conclude that \I/ 6 • KijS L s^ is not a good 
function to be prescribed in a CTS approach. The 
present numerical analysis suggests that adaptation of 
the coordinate system to the horizon (boundary condi- 
tion b = N) and the prescription of 9^ to a sufficiently 
negative value determines a unique solution. This is a 
relevant information for the dynamical 7 case, in particu- 



In order to implement quasi-equilibrium, and according to Eq. 
1141 . we restrain here to a flat 7jj (this guarantees the umbilical 

condition — H = 0). 7 In the dynamical case, the boundary condition b = N is gener- 
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Kijs's-i = A, k = 0.3 

-2.6 I 1 1 1 1 1 1 



, , -2.65 - 




-1.4 -1.3 -1.2 -1.1 -1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 



FIG. 5: Values of * 6 • K ij s i s i /H and 6 {k) /H, for 
boundary conditions KijS t s J — A, k — 0.3 and f2 = 
0.00, 0.02, 0.04, 0.18, 0.20. Curves with Q a / do not pass 
through the origin: conclusions in Ref. JT^| do not extend to 
the rotating case. On the other hand, £Iq 7^ curves become 
indistinguishable of the spherically symmetric curve (dotted 
line) for sufficiently negative values of KijS l s : > . 



lar in a constrained evolution scheme in which the elliptic 
system (|12H is solved at each time step (see Ref. |9()- 



C. Prescription of the divergence of flj 

Since the solution lnJVo of Eq. © is defined up to a 
constant on S, the value N = C ■ e lnN ° provides a Dirich- 
let condition for the lapse, which takes into account the 
manner in which S is embedded in the spacetime and per- 
mits to control the magnitude of N via the free choice 
of the constant C . The multiplicative character of the 
latter is a good feature on numerical grounds. Of course, 
condition (JSJ must be completed with a choice of 2 D k Ctk 
(an arbitrary function on S ~ S 2 with vanishing £ = 
mode). This makes this condition an effective one, in 
the sense discussed in Ref. 4] for "lapse boundary con- 
ditions". However, for specific problems natural choices 
exist, e.g. 2 D k ftk = for perturbations of the spherically 
symmetric case (see 0, El f° r a Pawlowski gauge, moti- 
vated by Kerr-Schild-like slicings). More importantly, a 
part of the extrinsic geometry is incorporated into this 
boundary condition, which thus goes beyond a merely 
numerically convenient choice. Once 2 D k Vlk{— 0) has 
been chosen, the slicing is completely fixed, and all the 
freedom is reduced to a single constant C that controls 
the rate to move across the slicing. We have numerically 
verified the existence of unique solutions when combin- 
ing this condition with those ones in the third column of 



Table |H 



D. General Comments 

Other possible combinations. Table[I]presents an associa- 
tion between boundary conditions and constrained fields. 
Even though this can be useful to organize the discussion, 
it must again be emphasized that this is only a method- 
ological option. Insisting on such an association can be 
misleading and can shadow some useful choices. As an 
example, since condition provides an expression for 
KijS l s^ , it could be interpreted as a condition 8 for b via 
Eq. I |lti[l . Then, interpreting b — N = as a condition 
for N, we can think of prescribing a value for \P on S by 
a Dirichlct condition. This particular example has two 
relevant applications: i) it permits to prescribe the area 
of S (since the conformal part is a free data) and ii) it 
provides an alternative strategy for the analytical study 
of the well-posedness of the marginally trapped condi- 
tion (Rcfs. |22l |23| focus on the control of the radial 
derivative of \& in order to apply a maximum principle to 
the Hamiltonian constraint equation; the alternative of a 
Dirichlct condition for 'F shifts the focus to the momen- 
tum constraint). 

Effective boundary conditions. At the end of the day, 
quasi-equilibrium conditions J3J together with b — N = 
leave a single function to be specified on S. Following 
El- IHl this fifth function can be seen as a condition 
"for the lapse" , either geometrically motivated or purely 
effective. But it must be underlined (see also in this sense 
the discussion in section III.C of Ref. 3), that it can 
be more generally interpreted as an effective condition 
on any scalar combination of fields, e.g. the conformal 
factor or 2 D k flk- 

Generic well-posedness. Finally, we have commented 
that combining conditions in Table [I] lead to well-posed 
problems for generic choices of f\ and fi. By this we 
mean that only for some critical f\ and fi the problem 
admits an infinite number of solutions. For instance, con- 
dition Ced(k) — in |2| can be recasted as a prescription 
fi for ^> 6 • KijS l sf and also as 2 D k flk = fz- If these par- 
ticular expressions are used for conditions (|16fl or (JSJl in 
the spherically symmetric case, this leads to an ill-posed 
problem, as shown in Ref. 0- By generic me mean here 
that a small perturbation of this critical cases makes the 
problem well-posed. 



alized by solving an elliptic equation on <S; see I3lll32l . 



In particular, together with 1141 as a condition on V 1 , this means 
that non-expanding conditions can be fully fulfilled by an appro- 
priate choice of the shift. 
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V. CONCLUSIONS 

This work represents the numerical counterpart of 
Refs. [Hl[l3)[l3> wnere isolated horizon boundary con- 
ditions were proposed. As a first result, the prescription 
of a constant value k — k on a instantaneously non- 
expanding horizon, using a coordinate system adapted 
to the horizon [i.e. conditions ©, © and Q], defines a 
well-posed problem. If k is set to K Kett {a, J), as proposed 
in the problem is degenerated (infinite number of 
solutions) in the spherically symmetric case and admits 
no solution when rotation is introduced. The only free- 
dom in this system is the choice of the constant n . The 
k = const condition does not enforce a quasi-equilibrium 
restriction on the geometry of the horizon. This means 
that it does not increase the physical degree of station- 
arity of the associated initial data (consistently with £| , 
where it is shown that physical quantities do not depend 
on the chosen lapse boundary condition). Its interest is 
rather in the evolution of such initial data, since it pro- 
vides a slicing where the lapse function must be initially 
time independent, something desirable numerically. 

Secondly, the results in Ref. ^| on the prescription of 
^P 6 • KijS l s : ' in the CTT system, do not extend straight- 
forwardly to the CTS case, except in spherical symme- 
try. In this particular case, there exists a negative bound 
Sasym such that for Sasym < ^ 6 ■ Kij s l / H < there is a 
unique solution. There exists a second bound S rn i n such 
that the prescription S m i n < ^S 6 ■ Kijs 1 s J ' / H < 5 asym ad- 
mits two degenerate solutions. For ty 6 ■ Kij s l / H < 5 m ; m 
no solution exists. In the rotating case, the strong depen- 
dence on the fifth boundary condition prevents us from 
deriving general bounds for \1/ 6 • KijS % sK However, for 
the particular choice of a coordinate system adapted to 
the horizon (b = N condition), the prescription of 6,^ 
to a sufficiently negative value guarantees the existence 
and uniqueness of a solution to the CTS elliptic system. 
Therefore, rather than ^ 6 • KijS l s J , the good parameter 
to be prescribed is Ofts . This represents an important in- 
formation for the implementation of evolving black holes 
as regular future trapping/dynamical horizons 0, @> 
in a constrained evolution scheme. 

Finally we have underlined the fact that prescrib- 
ing non-expanding and adapted-coordinate-system con- 
ditions leave one free function to be specified on the hori- 
zon. Due to the (non-linear) coupled nature of boundary 
conditions, this fifth condition is not specifically related 
to a particular field and, even though it can be useful to 
interpret it as a lapse boundary condition, other choices 
(e.g. a Dirichlct condition on \&) can prove to be useful. 
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APPENDIX A: COMMENTS ON FIGURE [2] 

In this appendix we derive some analytical relations 
that are valid for CTS data in the conformally flat, max- 
imal, spherically symmetric case. In particular we prove 
that the graph displayed in Fig. [5] is independent of a 
specific choice of inner boundary conditions on N and 
b (the prescription KijS^s 3 = X determine one point in 
the curve), and give a parametric analytical expression 
through which this curve is defined. 

a. Independency of Fig\^ . In spherical symmetry, 
with (r, 9, (f>) being spherical coordinates in which the ap- 
parent horizon is located at r = rjj, the maximal slicing 
CTS equations l|12f> [with a = 0] reduce to: 



d 2 ^ 2 d^ 

dr 2 r dr 
d 2 



12N 1 



d_ 

dr 



dr 2 (^) + ^(^) 



dr 2 



(A*) 



2 d(NV) 7* E 



dr 



12N 



(r- x /r) 



d 
dr 

d_ 
dr 





[log (A*- 6 )] 

2 



(Al) 
= 



= 



and the apparent horizon boundary condition (note that 
[3 e = {3<t> = 0) is given by 



dr 



2r 



d 

6JV* 







(A2) 



With the introduction of the compactificd radial coordi- 
nate s = rii/r, we get 



B 



12A 2 



= 



B" - B' 



s 2 (A*)" - IB 



log [A*~ 
vp 5 



12 



12N 2 



and 



~2 



6A 



B' 



0. 



where ' = d/ds and (3 r 
solved explicitly: 



-S=l 

B/s. Equation 



B' = c ■ ,s 2 A<3>~ 



(A3) 
(A4) 

(A5) 

(A6) 
can be 

(A7) 



where c is a constant of integration, closely related to 
* 6 • K lj s i s j [see Eq. JTJJ and below]. Using l|A7)l we 
obtain 



12 







(A8) 
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together with the boundary condition 



1 



c 

6*3 



(A9) 



s=l 



The system (|A8I IA9|l uniquely defines a sequence of solu- 
tions \I/(s;c) which is independent of a specific choice of 
inner boundary conditions for B and N . Moreover, the 
horizon quantity (cf. Eq. (|16[) m the spherically sym- 
metric case) 



sV] 

2 
3N 



r=r H 
d 

dr 



(A10) 



I3 r - -p 1 



where 



^ H = *(s= l;c) 



(All) 



is a function of c alone. Thus also the graph displayed 
in Fig. [3 is independent of a specific choice of inner 
boundary conditions on (3 r and N. The prescription 
KijS t s J = A determines a point in this curve. 

b. Analytic representation of Fig. In order to ob- 
tain explicit mathematical expressions that describe the 
graph displayed in Fig. [3 we consider the spatial metric 
as well as the extrinsic curvature of the family of time- 
independent maximal slicings of the Schwarzschild solu- 
tion 301, i.e. 



,/.*-= j I -^- + C 2 ^-) (/./)'"' -f U'-(tlH- + sin 2 



Rh 



and 



R2 

k 1 - r -H- 



which gives in particular 




2C 
Rh 



(A12) 
(A13) 

(A14) 



In these coordinates (R, 6, <fi), the radius of the black hole 
horizon is given by R = Rh = 2M, where M is the black 
hole mass. The constant C parametrizes the family of 
maximal slicings. 

The coordinate transformation r — r(R), leading to a 
conformally flat line element 

ds 2 = * 4 [dr 2 + r 2 (d0 2 + sin 2 9d(j) 2 )] , (A15) 

is described by 



s = exp 



da 



/i ay/l - a + C 2 a 4 / ' 



(A16) 



where 



s = — and S = — ^- 
r R 



The conformal factor 4" which satisfies the boundary 
value problem QA8I IA9|I . can be obtained from the com- 
parison of the line elements l)A12l IA15|) : 



> R Rh s 
r rn S 



As a consequence 



y H — 



Rh 



(A18) 



(A19) 



and by imposing the asymptotic boundary condition, 
lims^o * = 1, we can write ^ 2 H in terms of the parame- 
ter C: 



1 



* 2 K (C>0)=e X p(7 ^ 



y/l-a + C 2 cr 4 



(A20) 

The formula l|A20(l only describes ^fjj for non-negative 
values of C. An expansion of this expression into the 
realm of negative C-values yields, for > C > — 3\/3/16 



* 2 H (C<0) = cxp^ 1 ^ 



da 
c 

si(C) 



1 



VI - (r + C 2 a 4 
da \ 



(A21) 



aVl-cr + C 2 a 4 r 
where si(C) is the real zero of 

f(s) = l- S + C 2 s A , 

with 

l<*i(C)<|. 

Note that ^ 2 H {C) tends to +oo as C -> -3V3/16 since 
f(s) has a double zero, s\ — 4/3, in this limit. 

With l|A14|) and H — 2/th, we finally obtain the de- 
sired parametric description of the curve displayed in Fig. 
121 [a parametrization in terms of c in Eq. (|A11I) can be 
also obtained from the relation between C and c provided 
by Eqs. (|A"T0|) and (£Q}]: 



Abscissa: K ljS l s^/H = -C^ H 2 
Ordinate: ^> 6 H ■ K ij s i s j /H = 



(A22) 



A particular consequence of this analysis is 



lim V b H - i^-sV 



-2H. 



(A23) 



(A17) providing the asymptotic value when KijS l s 3 / H — > — oo. 
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